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The effect of power interchange operation of multiple household gas engine cogeneration systems (H- 
GCGS) on the energy-saving is investigated using an optimization approach based on the mixed-integer 
linear programming. In this power interchange operation, electricity generated by H-GCGS is shared 
among households in a housing complex without transmitting to a commercial electric power system so 
that the operating time of these systems may increase. This paper numerically analyzes optimal oper¬ 
ational strategies for 20 households and three types of household energy supply configurations: the 
power interchange operation of the H-GCGSs (IC), stand-alone operation of each H-GCGSs (SA), and 
conventional energy supply system without the H-GCGSs. A numerical result clarifies the effectiveness of 
the power interchange operation from the energy-saving viewpoint and a dominant parameter for 
evaluating the energy-saving effect. 

© 2008 Elsevier Ltd. All rights reserved. 


1. Introduction 

Cogeneration systems have been increasingly utilized for 
industrial and commercial applications because of their high 
energy-saving and economic efficiencies. Owing to recent devel¬ 
opments of small-scale prime movers with high performance 
including gas engines (GE), fuel cells, and Stirling engines, their 
applications have been extended to households [1]. Most of oper¬ 
ating household cogeneration systems adopt GEs because of their 
low cost and high reliability. 

The energy-saving effect by utilizing household gas engine 
cogeneration system (H-GCGS) has been studied. Paepe et al. 
reported the energy-saving effect of 4.7 and 5.5 kW H-GCGSs 
that are operated at a detached house and terraced family house 
in Belgium [2]. Cockroft and Kelly reported that an H-GCGS 
meeting the maximum heat demand in a household needs to 
operate with overall efficiency of more than 80% for energy¬ 
saving [3]. Furthermore, Hawkes and Leach reported the C0 2 
reduction effect through an analysis of optimal operational 
strategy of 2 kW H-GCGS [4]. In these studies, the H-GCGSs can 
transmit surplus electricity to commercial electric power 
systems; thus, the H-GCGSs can operate to meet heat demand 
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variations and do not need storage tanks (STs). However, in 
Japan, electricity generated by household cogeneration systems 
cannot be transmitted to commercial electric power systems by 
the regulation. Therefore, H-GCGSs for Japanese household are 
downsized and generate electricity at a rating of approximately 
1 kW [5,6]. These systems have an operational restriction that 
they can operate only when the electricity demand of each 
household exceeds the generated electricity at the rating. This is 
because they must adopt on-off operation that implies either 
rated load operation or shutdown due to the inefficiency of their 
partial load operations [7]. Moreover, a ST must be installed for 
intermittent hot water supply under this operational restriction; 
however, heat loss from the ST is considerable when the oper¬ 
ating period of the H-GCGS is much different from the high 
demand period of hot water. Therefore, actual energy-saving 
effect of the H-GCGSs in Japan is considered to underachieve. To 
overcome this problem, the capacity of the GE can be further 
reduced; however, this is not effective because generating effi¬ 
ciency of a smaller scale GE markedly decreases. 

Based on this background, the authors propose power inter¬ 
change operation of multiple H-GCGSs in a housing complex. In this 
operation, electricity generated by the H-GCGSs is shared among 
households in the housing complex without transmitting to an 
external commercial electric power system. Since the surplus 
electricity generated at a household with little electricity demand 
can be transmitted to other households that have much electricity 
demand, the aforementioned their operational restriction is 
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Nomenclature 

E 

electricity [kWh/h] 

F 

natural gas consumption [m 3 /h] 

j 

objective function (daily primary energy 
consumption) [MJ/d] 

I< 

number of sampling times 

k 

index for sampling times 

M 

number of representative days 

m 

index for representative days 

N 

number of object households 

n 

index for object households 

p, q, r, s 

performance characteristic values [l<Wh/h,m 3 /h, 
(MJ/h)/(m 3 /h), kWh/h] 

a 

heat flow rate [kWh/h] 

S 

stored heat [MJ] 

At 

sampling time interval [h] 

a 

reduction rate in total daily primary energy 
consumption [%] 

5 

binary variable for expressing on-off status 

Vi 

heat loss rate of storage tank [%/h] 

0 

upper limit 

Abbreviation (subscript) 

CO 

conventional energy supply system 

GE 

gas engine 

GW 

gas-fired water heater 

H-GCGS 

household gas engine cogeneration system 

IC 

power interchange operation of household gas 
engine cogeneration systems 

SA 

stand-alone operation of each household gas engine 
cogeneration system 

ST 

storage tank 

THPR 

total daily heat to power demand ratio of object 
households 

Superscript 

a 

auxiliary machine 

in 

inlet 

out 

outlet 


relaxed. As a result, each H-GCGS under the power interchange 
operation can operate to meet heat demand variations and to 
minimize the heat loss from the ST. In addition, there is hardly any 
electric transmission loss in the housing complex. Hence, the 


power interchange operation is expected to increase their energy¬ 
saving effect. 

This paper discusses the energy-saving effect of the power 
interchange operation as the first step of its feasibility study. To 
clarify theoretical limit of the energy-saving effect, operational 
strategies of multiple H-GCGSs are numerically analyzed by an 
optimization approach based on the mixed-integer linear 
programming. The analysis is conducted for 20 households and 
three types of household energy supply configurations: the afore¬ 
mentioned power interchange operation of the H-GCGSs, stand¬ 
alone operation of each H-GCGS, and conventional energy supply 
system without the H-GCGS. Through a numerical analysis of the 
operational strategies, effectiveness of the power interchange 
operation from the energy-saving viewpoint, and a dominant 
parameter for evaluating the energy-saving effect are clarified. 

2. System configuration 

Fig. 1 shows a schematic diagram of power interchange opera¬ 
tion of multiple H-GCGSs in a housing complex (IC). Each household 
in a housing complex installs the H-GCGS, a ST, and a gas-fired water 
heater (GW). The solid, dash-dotted, and dotted lines in this figure 
show the flows of hot water, electricity, and natural gas, respec¬ 
tively. This study considers electricity and hot water as energy 
demands of the households because space cooling and heating are 
supplied by electric air conditioners installed at each household. 

In IC, electricity generated by the H-GCGSs is shared among the 
households in the housing complex; however, the generated elec¬ 
tricity is not transmitted to an external commercial electric power 
system. Thus, electricity is supplied to meet the electric demands of 
the households by operating the H-GCGSs and by purchasing 
electricity in bulk from an electric power company. Generated 
electricity is also provided to drive the auxiliary machines of the H- 
GCGSs. This power interchange operation has little influence on the 
stability of the electric power system because electricity generated 
by each H-GCGS is synchronized with the electric power system 
through an inverter for voltage and phase conditioning. Exhaust 
heat from both jacket cooling water and exhaust gas of the GE is 
recovered as hot water and stored in the ST. Hot water in the ST is 
independently consumed at each household. The shortage of hot 
water from the ST is supplemented by the GW. The H-GCGSs are not 
operated with a heat disposal from the energy-saving viewpoint. In 
actual operation, a control system, which manages operation of 
multiple H-GCGSs in response to variations in electricity and hot 
water demands of object households and stored heat in STs, and 
on-off accessory for each H-GCGS are required, although they are 


Commercial 
electric power 
system 



Fig. 1 . Schematic diagram of power interchange operation of H-GCGSs. 
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not considered in this study because the energy-saving effect is 
evaluated for the feasibility study of IC. Thus, on-line control 
method for IC must be developed separately for the practical use. 

This study also considers two other types of household energy 
supply configurations to evaluate the energy-saving effect of IC. 
One is stand-alone operation of each H-GCGS without the power 
interchange among the household (SA). This is the only operation 
method of the H-GCGS permitted in Japan. The other is conven¬ 
tional energy supply system without the H-GCGS (CO): electricity is 
supplied only by purchased one from the electric power company, 
and hot water is supplied only by the GW. 

3. Optimal operational planning problem 


denote the natural gas consumption of the H-GCGS and GW, 
respectively; and Ep denotes the electricity purchased in bulk. For 
SA, the objective function Jsa is expressed by the following 
equation: 


JsA(m) 


I< 

ftE 


k = 1 


N N 

J2 f ge(", k,m) + J2 F cw (n, k, m) 


n = 1 


n = 1 


At 


I< N 

+ <Pe EE Ep(n,/c, m)At (m 

/<= 1 n = 1 




where E\> denotes the electricity purchased individually at each 
household. For CO, the objective function Jco is expressed by the 
following equation: 


The energy-saving effect of IC is evaluated using an optimization 
approach based on the mixed-integer linear programming. That is, 
the optimal operational strategies of the three types of household 
energy supply configurations including on-off and rated-partial 
load status of the operation of the system components, energy flow 
rates, and stored heat are determined to minimize primary energy 
consumption subject to the satisfaction of energy demand 
requirements. This optimal operation planning problem is 
described below. 


3.1. Mathematical formulation 

The housing complex investigated in this study has N house¬ 
holds, and the index for the households is designated by n, i.e., 
n = 1,2,..., N. Moreover, to consider seasonal and hourly changes in 
the energy demands, a typical year is divided into M representative 
days. The index for the representative days is designated by m, i.e., 
m = 1,2,..., M. Furthermore, each representative day is divided into 
I< sampling times with a period of At, i.e., At = 24//C The index for 
the sampling times is designated by k, i.e., k = 1, 2, ..., I<. The 
operational strategy of a cogeneration system with a ST cannot be 
determined independently for each sampling time. Thus, the daily 
cyclic operation, which assumes that energy demands change 
cyclically with a period of 24 h on each representative day, is 
considered [8]. 


3.1.1. Decision variables 

The decision variables are composed of binary and continuous 
variables. The binary variable expresses the on-off status of the 
operation of the H-GCGS of the nth household at the kth sampling 
time on the mth representative day. The continuous variables 
express the flow rates of the input and output energies of the 
system components and stored heat in the ST in the same condition 
as that of the binary variables. 


3.1.2. Objective function 

The objective function J to be minimized is the total daily 
primary energy consumption of the object households on each 
representative day in order to evaluate the energy-saving effect. 
The definition of the total daily primary energy consumption 
depends on the household energy supply configurations. 

First, for IC, the objective function Jic on the mth representative 
day is expressed by the following equation: 


K 


Jic(m) = <p c E 
k= 1 


N N 

E f CE(n, km)+J2 F cw(n, k, m) 

n = 1 n = 1 


At 


I< 


+ <PE E*p b (/<, m)At (m = 1,2, ...,M) 
k=\ 


( 1 ) 


where cpc and cp E denote the conversion factors for the primary 
energy of natural gas and electricity, respectively; Fge and Few 


I< N 

Jco(m) = ftE E F cw (n,k,m)At 

k = 1 n = 1 
K N 

+ , e E E4(n,k,m)At (m = 1,2, (3) 

k = 1 n = 1 


3.1.3. Constraints 

The constraints in this optimization problem are composed of 
the performance characteristics of the system components and 
energy balance and supply-demand relationships. 

The H-GCGSs adopt the on-off operation, which implies either 
rated load operation or shutdown [2,7]; this is attributed to the 
characteristic that the partial load operations of the small-scale GEs 
are inefficient. Thus, the generated electricity Fge, heat flow rate of 
recovered hot water Qge, natural gas consumption Fge, and elec¬ 
tricity consumed by the auxiliary machine F^ E are expressed as 
follows: 


E GE (n,k,m) = p GE d CE (n,k,m) ' 
0c E (n,k,m) = q CE 5 GE {n,k,m ) 
F GE (n,k,m) = r GE d GE (n,k,m) > 
E GE (n,k,m) = s GE d CE (n,k,m ) 
<5 GE (n,k,m)e{0,l} 



12 N- 


k = 1,2, m 




where <Sge denotes the binary variable expressing the on-off status 
of the operation of the H-GCGS; the coefficients Pge, 9ge, Jge, and Sge 
denote the performance characteristic values of the H-GCGS. 

For the STs, the heat balance is considered by the following 
difference equation: 


S ST (n,/<,m)-S ST (n,/<-l,m) ^ in/ . . 

^ v — L —f = Q^(n,/c, m) 

-71 L AtS ST (n,k- l,m) 

(n = 1,2, ...,7V; k = 1,2, ...,K; 
m = 1,2, ...,M) 


Qsf lt (n 7 k, m) 



where Sst denotes the heat stored in the ST; and denote 
the heat flow rates of hot water stored into and supplied from 
the ST, respectively; 77 l denotes the heat loss rate of the ST. 
Moreover, the heat stored in the ST is limited by its capacity S ST 
as follows: 


0 < S ST (n,k,m) < S ST (n = 1,2, ...,N; k = 0,1, 

m = 1,2,...,M) (6) 

Under the assumption that the ST is cyclically used with 
a period of 24 h, the heat stored at the initial state is equal to 
that at the terminal state at each household on each represen¬ 
tative day. 
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Table 1 

Performance characteristic values of system components 


Equipment 

Performance 
characteristic value 

Gas engine cogeneration system 

Pge =1.0 kW 


Qge = 3.7 kW 


r GE = 0.443 m 3 /h 


s GE = 0.075 kW 

ST 

S ST = 30.1 MJ 


ijl= 1.0%/h 

GW 

<Jcw = 36.0 MJ/(m 3 /h) 


S ST (n,0,m) = S ST (n,I<,m ) (n = 1,2 m = 1,2, 

(7) 

For the GWs, the relationship between the heat flow rates of 
supplied hot water Qgw and the natural gas consumption Few are 
expressed by the following equations: 

Qcw(n, k, m) = q G w F cw (n, k,m)\ 

k = 1,2m = 1,2(8) 

where the coefficient gew denotes the performance characteristic 
value of the GW. 

Furthermore, the energy balance and supply-demand rela¬ 
tionships for electricity and hot water are considered. The 


electricity balance and supply-demand relationships depend on the 
household energy supply configurations: for IC, it is considered in 
the housing complex; for SA and CO, it is considered at each 
household. On the other hand, the hot water balance and supply- 
demand relationships are considered at each household in any 
household energy supply configuration. 


3.2. Solution method and its validity 

The optimal operational planning problem formulated in the 
previous section results in a mixed-integer linear programming 
problem. This is solved by the GAMS (General Algebraic Modeling 
System)/CPLEX solver, which is a general-purpose solver that 
combines the branch and bound method with the simplex one for 
large-sized mixed-integer linear programming problems [9]; thus, 
optimal operational strategies calculated by this solver have high 
accuracy. 


4. Numerical study 

A numerical study on the energy-saving effect of IC is conducted 
using the aforementioned optimization model. First, the input data 
for this optimization problem are described; then, the energy¬ 
saving effect of IC is discussed through a numerical analysis of 
optimal operational strategies. 
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Fig. 2. Hourly electricity demands at 20 households in December, (a) House-1 to house-10, (b) house-11 to house-20. 
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Fig. 3. Hourly hot water demands at 20 households in December, (a) House-1 to house-10, (b) house-11 to house-20. 


4.1. Input data 

The input data for this optimization problem are composed of 
the performance characteristic values of the system components, 
household energy demands, and conversion factors for primary 
energy. 

The performance characteristic values of the system compo¬ 
nents are listed in Table 1. These values are derived from Refs. [7] 
and [10] and our actual operating performance. The net generating 
efficiency considering the electricity consumed by the auxiliary 
machine and heat recovery efficiency of the H-GCGS and thermal 
efficiency of the GW calculated using the higher heating value of 
natural gas are about 17%, 67%, and 80%, respectively. 

The housing complex investigated in this study has 20 house¬ 
holds, i.e., N = 20. The energy demands at the 20 households are 
estimated at 24 sampling times for a representative day of each 
month, i.e., M=12, K= 24, and At=l h. The electricity and hot 
water demands at the 20 households have been measured by Tsuji 
et al. [ 11 ]. As examples of the energy demand patterns, Figs. 2 and 3 
show the hourly electricity and hot water demands at the 20 
households on the representative day of December, respectively. 


Table 2 

Conversion factors for primary energy 


Energy source 

Value 

Electricity 

q> E = 9.83 MJ/kWh 

Natural gas 

cpc = 45.0 MJ/m 3 


The hourly electricity demands differ widely among the 20 
households. Moreover, the period when the electricity demand 
exceeds the net generated electricity of the H-GCGS is limited at 
any household. On the other hand, the hot water demands are 
concentrated in the nighttime at any household. 

The conversion factors for primary energy are listed in Table 2. 
For electricity, the thermal power average conversion factor defined 
in the Law Concerning the Rational Use of Energy of Japan is 
introduced. 

4.2. Results and discussion 

4.2.1. Energy-saving effect of power interchange operation for 
five households 

The energy-saving effect of IC is discussed through the numer¬ 
ical analysis of the optimal operational strategies. First, the energy¬ 
saving effect of IC for five of 20 object households, which are called 
as House-1 to -5 for convenience, i.e., N = 5, is investigated. Their 
energy demands are shown in Figs. 2 and 3. 

Figs. 4 and 5 show the optimal operational strategies using IC 
and SA on the representative day of December, respectively. The 
optimal operational strategy for the electricity supply using IC 
shows the total electricity supply to the five households, and the 
other optimal operational strategies show the energy supply only 
to House-3. As shown in Fig. 5, the H-GCGS at House-3 cannot be 
operated in SA because the electricity demand at House-3 is less 
than the net generated electricity of the H-GCGS throughout the 
day. Therefore, all hot water is supplied by the GW. On the other 
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Fig. 4. Optimal operational strategy using IC in December, (a) Total electricity supply, 
(b) hot water supply (house-3). 


hand, the H-GCGS at House-3 in IC can operate for 4 h because the 
total electricity demand of the object households becomes larger 
than the net generated electricity of the H-GCGS. Thus, IC relaxes 
the operational restriction of the H-GCGS in SA and increases the 
operating time of the H-GCGSs. In this case, most of hot water is 
supplied from the ST. Moreover, the operation of the H-GCGSs in IC 
is concentrated in the nighttime so that the heat losses from the STs 
may be minimized. This operational strategy contributes to the 
reduction in the natural gas consumption for the GW. 

Table 3 lists the computational result of the total daily 
purchased electricity, and the total daily natural gas and primary 
energy consumptions of the five households on the representative 
day of December for the three types of household energy supply 
configurations. By the introduction of the H-GCGSs (IC, SA), the 
total daily natural gas consumption of the five households increases 
whereas the total daily purchased electricity decreases. These are 
because of operation of the H-GCGSs. As a result, the total daily 
primary energy consumptions in IC and SA are less than that in CO. 
Furthermore, the total daily primary energy consumption in IC is 
lower than that in SA owing to the increase in the operating time of 
the H-GCGSs in IC. 

4.2.2. Effect of increasing number of object households 

On the basis of the result for the five households, the energy¬ 
saving effect due to increasing the number of object households is 
investigated. 

Fig. 6 shows the optimal operational strategies for the total 
electricity supply using IC on the representative day of December at 
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Fig. 5. Optimal operational strategy using SA at House-3 in December, (a) Electricity 
supply, (b) hot water supply. 


10 (House-1 to -10), 15 (House-1 to -15), and 20 households 
(House-1 to -20), i.e., N= 10, 15, 20. For any number of the object 
households, the operation of the H-GCGSs is concentrated in the 
nighttime when the hot water demand increases. Moreover, by 
analyzing the optimal operational strategies for the hot water 
supply, it is confirmed that most of hot water is supplied from the 
ST at any household. 

To evaluate the energy-saving effect of IC, the reduction rates in 
the total daily primary energy consumption from CO to IC ctic/co and 
those from SA to IC chc/sa for four cases of the number of the object 
households, i.e., N= 5, 10, 15, 20, on the representative days of 
August, October, and December are listed in Table 4. ctic/co and cuc/sa 
on the mth representative day are defined as follows: 


a \c/co( m ) 


Jco(m) 


X 100 





a IC/SA( m ) 


JsA(m) -Jic(m) 

JsA(m) 


X 100 




( 10 ) 


Table 3 

Total daily purchased electricity, and total daily natural gas and primary energy 
consumptions of five households in December 


Energy supply configuration 

IC 

SA 

CO 

Total purchased electricity, kWh/d 

73.3 

81.6 

92.7 

Total natural gas consumption, m 3 /d 

9.87 

9.18 

8.23 

Total primary energy consumption, MJ/d 

1165 

1215 

1282 
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Fig. 6. Optimal operational strategy for total electricity supply using IC in December, (a) 10 households, (b) 15 households, (c) 20 households. 


where Jco, Jsa, and J\c denote the total daily primary energy 
consumption in CO, SA, and IC, respectively. 

Obviously, IC has a higher energy-saving effect than SA and CO on 
any representative day and for any number of the object households. 
The increase of ctic/co is outstanding in December where the oper¬ 
ating time of the H-GCGSs increases strikingly owing to more hot 
water demand, while ctic/sA increases in October and December. 
However, there is no definite correlation between the two types of 
the reduction rates and the number of the object households. The 


aforementioned result of the optimal operational strategy for the 
five households indicates that the energy-saving by IC mainly results 
from the increase in the operating time of the H-GCGSs at the 
households that have less electricity demand. Thus, the energy¬ 
saving effect of IC is not affected by the number of the object 
households but affected by the energy demand patterns of the object 
households. 

To clarify a dominant parameter to characterize the energy¬ 
saving effect of IC, the relationship between the two types of the 
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Table 4 

Reduction rate of total daily primary energy consumption for four cases of number of 
object households on three representative days 


Reduction rate of total daily primary energy consumption a % 

Month 

August 


October 


December 

Comparative configuration 

IC/CO 

IC/SA 

IC/CO 

IC/SA 

IC/CO 

IC/SA 

5 households 

3.32 

1.08 

6.65 

3.16 

9.14 

4.17 

10 households 

2.60 

1.02 

6.38 

3.60 

9.33 

3.51 

15 households 

2.70 

1.19 

6.48 

4.36 

9.64 

4.42 

20 households 

2.86 

1.22 

7.08 

5.52 

10.04 

5.62 


reduction rates in the primary energy consumption and the total 
daily heat to power demand ratio of the object households (THPR) 
is examined. THPR on the mth representative day is defined as the 
ratio of the total daily hot water demand of the object households 
to their total daily electricity demand as follows: 

K N 

E E Qo(n,k,m)At 

THPR(m)=^l^- (in = 1,2,...,M) (11) 

E E E D (n,k,m)At 

k= 1 n = l 

where £ D and Qd denote the demand of electricity and hot water of 
the nth household at the kth sampling time on the mth represen¬ 
tative day, respectively. 

The numerical result of the relationship based on the optimal 
operational strategy on the representative day of each month is 
shown in Fig. 7. a\qco increases almost linearly with THPR; there¬ 
fore, the energy-saving effect by introducing IC can be evaluated 
using THPR. On the other hand, a\qsA tends to increase with THPR; 
however, there is no clear correlation between cqc/sA and THPR. 
Since the relaxation of the operational restriction by IC is affected 
by the energy demand patterns of the object households, a detailed 
analysis of the energy demand patterns is required to investigate 
a dominant parameter to characterize cqc/sA that is the criterion to 
determine the energy supply configuration using the H-GCGS. 

4.2.3. Difference between actual energy-saving effect and 
optimal result 

The abovementioned optimal results are computed under the 
condition where the energy demands of the object households are 
preliminarily determined. This means that the predicted result of 
the household energy demands totally coincides with the actual one. 
Thus, the energy-saving effect of the H-GCGSs under actual opera¬ 
tion may be lower than the optimal result due to uncertainty of the 
energy demand prediction. However, it is considered that the 
decrease in the actual energy-saving effect in IC may be less than that 
in SA. This is because a prediction error rate of the total electricity 
demand of object households, which becomes the electric load to 
each H-GCGS in IC, can be smaller than that of each household owing 
to offsetting of prediction error at each household. 

Furthermore, the sampling period of the energy demands At 
also has a significant influence on the energy-saving effect of the 
H-GCGSs derived from an optimization model because actual 
household energy demands generally have fluctuations with short 
periods due to the use of appliances. Hawkes and Leach reported 
that CO 2 reduction effect of household cogeneration systems in SA 
is overestimated by analyses conducted using the sampling period 
of 1 h [12]. Thus, the energy-saving effect of SA in this study may be 
also overestimated; however, it is considered that the over¬ 
estimation of the energy-saving effect of IC may be smaller than 
that of SA. This is because the fluctuation of electricity demand at 
each household is smoothed by evaluating the total electricity 
demand of the object households in IC. 

Therefore, IC is expected to have an advantage in the decrease in 
the energy-saving effect under actual operation over SA. 
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Fig. 7. Relationship between reduction rate in total primary energy consumption and 
THPR in each month, (a) a IC /co, (b) «ic/sa- 


5. Conclusions 

The energy-saving effect of power interchange operation of 
multiple household gas engine cogeneration systems (H-GCGSs) in 
a housing complex has been investigated by an optimization 
approach as the first step of its feasibility study. In the power 
interchange operation, electricity generated by the H-GCGS is 
shared among households in the housing complex without trans¬ 
mitting to an external commercial electric power system. A 
numerical analysis of optimal operational strategies by the mixed- 
integer linear programming has been conducted for 20 households 
and three types of household energy supply configurations: the 
power interchange operation of the multiple H-GCGSs (IC), stand¬ 
alone operation of each H-GCGSs (SA), and conventional energy 
supply system without the H-GCGSs (CO). The results obtained in 
this study are summarized as follows: 

• IC increases the operating time of the H-GCGS and enables their 
operation so that the heat loss from the ST may be minimized. 

• IC is the most effective in the energy-saving among the three 
types of the household energy supply configurations in any 
season for any number of the object household. 

• There is no correlation between the energy-saving effect of IC 
and the number of the object households. 

• The energy-saving effect from CO to IC can be evaluated with 
the THPR. 
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These results clarify the energy-saving effect of IC in the housing 
complex in the area where electricity generated by household 
cogeneration systems cannot be transmitted to commercial electric 
power systems. As future tasks, the authors will investigate 
a dominant parameter to characterize the energy-saving effect of IC 
against SA through a detailed analysis of the energy demand 
patterns of the object households. Furthermore, the authors 
recognize the importance of financial benefit by the H-GCGSs. 
However, in order to conduct the feasibility study from the 
economic viewpoint, the information including the periodic 
maintenance cost, additional cost of a control system for the power 
interchange operation, and utility charge structure for the power 
interchange operation in a housing complex as well as the initial 
cost of the H-GCGS must be required. Thus, this feasibility study is 
considered as the second step of our study. 

References 

[1 ] D’Accadia MD, Sasso M, Sibilio S, Vanoli L. Micro-combined heat and power in 
residential and light commercial applications. Appl Therm Eng 
2003 ;23( 10): 1247-59. 

[2] Paepe MD, D’Herdt P, Mertens D. Micro-CHP systems for residential applica¬ 
tions. Energy Convers Manage 2006;47(18-19):3435-46. 


[3] Cockroft J, Kelly N. A comparative assessment of future heat and power 
sources for the UK domestic sector. Energy Convers Manage 2006;47(15- 
16):2349-60. 

[4] Hawkes AD, Leach MA. Cost-effective operating strategy for residential micro- 
combined heat and power. Energy 2007;32(5):711-23. 

[5] Osaka Gas Co., Ltd. Annual report 2003. Osaka, Japan: Osaka Gas Co., Ltd.. 
Available from: <http://www.c-direct.ne.jp/osaka/english/contents/irlibrary/ 
pdf/ar2003e/ar2003e.pdf>; 2003. 

[6] Osaka Gas Co. Ltd, Toho Gas Co. Ltd., Saibu Gas Co. Ltd., Honda Motor Co. Ltd., 
Noritz Co. Ltd., Chofu Seisakusho Co. Ltd. Development of “ECOWILL” cogen¬ 
eration system for residential applications. Japan gas association newsletter. 
Tokyo. Japan: Japan Gas Association. Available from: <http://www.gas.or.jp/ 
english/pdf/ecowill.pdf>; 2004. 

[7] Osaka Gas Co. Ltd. Press release (y° ^ IJ I ) —^Zj. Osaka, Japan: Osaka Gas 

Co. Ltd.. Available from: <http://www.osakagas.co.jp/Press/pr02/021226. 
htm>; 2002 [in Japanese]. 

[8] Yokoyama R, Ito K. Optimal operational planning of cogeneration systems with 
thermal storage by the decomposition method. ASME J Energy Resour Tech 
1995;117(4):337-42. 

[9] Brooke A, Kendrick D, Meeraus A, Raman R. GAMS a user’s guide. Washington, 
DC: GAMS Development Corp.; 1999. 

[10] Honda Motor Co. Ltd. Press release. Tokyo, Japan: Honda Motor Co. Ltd.. 
Available from: <http://world.honda.com/news/2004/c041020.html>; 2004. 

[11 ] Tsuji Laboratory. Monitored energy data for residential houses in kansai region 
1998-2005. Osaka, Japan: Graduate School of Engineering Osaka University, 
Division of electrical, electric and information engineering: 2007. 

[12] Hawkes AD, Leach MA. Impacts of temporal precision in optimisation 
modelling of micro-combined heat and power. Energy 2005;30( 10): 1759-79. 


